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RESEARCH MEMORANDUM 

A CORRELATION BY MEANS OF THE TRANSONIC SIMILARITY RULES 
OF THE EXPERIMENTALLY DETERMINED CHARACTERISTICS OF 
22 RECTANGULAR WINGS OF SYMMETRICAL PROFILE 
By John B. McDevitt 

SUMMARY 


The transonic similarity rules have been applied to the correlation 
of experimental data for a series of 22 rectangular wings having sym- 
metrical NACA 63A-series sections, aspect ratios from l/2 to 6, and 
thicknesses from 2 to 10 percent- The data were obtained by use of the 
transonic bump technique over a Mach number range from 0.40 to 1.10, 
corresponding to a Reynolds number range from 1.25 to 2.05 million. 

The results show that it is possible to correlate experimental data 
throughout the subsonic, transonic, and moderate supersonic regimes by 
using. the transonic similarity parameters in forms which are consistent 
with the Prandtl-Glauert rule of linearized theory. ’ 

The multiple families of basic data curves for the various aspect 
ratios and thickness ratios have been summarized in single presentations 
involving only one geometric variable - the product of the aspect ratio 
and the l/3 power of the thickness ratio. 


INTRODUCTION 

A unified approach toward an understanding of transonic flows has 
been achieved only in recent years. Our meager knowledge of transonic 
flows, in comparison ■with the more complete and cogent understanding of 
subsonic and supersonic flows, is due not only to the complexities of 
the mathematics involved but also to the limitations of test facilities 
at transonic speeds. 

Similarity rules for transonic flow in two dimensions were derived 
by von Kerman (reference 1) and were extended recently by Spreiter 
(reference 2) and Berndt (reference 3) to include wings of finite span. 
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Biese r ul es can "be shown bo coincide with, one of the possible forms of 
the Prandtl-Glauert rule of subsonic and supersonic flows. (Reference 2.) 
Although the similarity rules do not provide explicit solutions, they do 
suggest the manner in which experimental data can "be correlated. 

A recent systematic experimental investigation of the effects of 
wing aspect ratio and thickness at transonic speeds (reference It-) has 
provided experimental data ideally suited to correlation using the tran- 
sonic similarity parameters. The analysis of these data is presented in 
this paper to provide the transonic characteristics of rectangular wings 
of symmetrical profile and to help evaluate the usefulness of the tran- 
sonic similarity parameters for the correlation of experimental data. 

The similarity rules are presented in slightly modified forms to 
permit a direct and convenient application in the data, correlation. The 
form of the data presentation waB chosen so that direct comparisons with 
the various linearized theories could be indicated. 

The correlation is applied to the experimental data for 22 rectan- 
gular wings having symmetrical MCA 63A -series sections, aspect ratios 
from l/2 to 6, and thicknesses from 2 to 10 percent. She data were 
obtained by use of the transonic bump technique over a Mach number range 
of 0.40 to 1.10, corresponding to a Reynolds number range from 1.25 to 
2.05 million. The basic data for the aspect ratios 1, 2, t, and 6 have 
been published in reference k. The basic data for the aspect ratios 0.5, 
1.5, and 3 were obtained by identical testing procedures but have not 
been published previously. A description of the wing models and the 
testing procedure is given in. reference 4. 

SYMBOLS' 


aspect ratio 
wing span 
wing chord 
total drag coefficient 
minimum drag coefficient 
minimum pressure-drag coefficient 

friction-drag coefficient, assumed equal to minimum drag 
coefficient at 0.7 Mach number 


A 

b 


CD 

c Dmin 

(CD P>min 

CD f 
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ac d 

drag coefficient due to lift (Cp - Cp^^) 

cl 

lift coefficient ( ) 

V as/ 

CLopt 

lift coefficient for maxi mum lift-drag ratio 

Cm 

pitching -moment coefficient, referred to 0 . 25 c 
f pitching moment^ 

\ <lSc / 

C p 

pressure coefficient 
ordinate distribution function 

<!■$ 

thickness distribution function 

<H) 

camber distribution function 

h 

camber parameter proportional to the amount of camber 

(d)^ 

maximum lift -drag ratio 

M 

free -stream Mach number 5 

Mcr 

critical Mach number 

Mpp 

drag-divergence Mach number 

<1 

dynamic pressure 

s 

wing area 

t/e 

thickness -to-chord ratio 

x,y,z 

Cartesian coordinates where x extends in the direction of 
the free-stream velocity 

a 

angle of attack 

V 

ratio of specific heats (for air y. = 1.10 

T 

ordinate -amplitude parameter 


k 
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<p perturbation-velocity potential 

cLCt 

— h , Cj q slope of lift curve, measured at zero lift coefficient 


slope of pitching-moment curve 

pressure -drag function 
- drag -due -to -lift function 
lift function 
pitching-moment function 
pressure-coefficient function 


THEORETICAL CONSIDERATIONS 
Transonic Similarity Rules 

Similarity rules (references 1, 2, and. 3) for the transonic flow 
around thin wings show how a wing can be affinely distorted in order to 
retain a related potential distribution when the free-stream Mach number 
is altered. These rules of correspondence can be shown to coincide with 
one of the possible forms of the ' Prandtl-Glauert rule, derived from the 
linearized potential equation of subsonic or supersonic flows (refer- 
ence 2). The various results of linearized theory can be expressed in 
terms of the transonic similarity parameters and the use of these param- 
eters will be consistent for both the linearized theory and the nonlinear 
transonic theory. 

Experience has shown that the linearized potential equation is a 
good approximation for sufficiently thin wings when the whole field of 
flow is subsonic or well-established supersonically. In the transonic 
regime certain nonlinear terms must be retained to show the mixed nature 
of the flow. The nonlinear terms of the basic potential equation are 
Important only near sonic speed and only the nonlinear terms involving 
derivatives with respect to the direction of the undisturbed fluid motion 
need be retained, since perturbations in sonic flow diminish more slowly 
in the directions perpendicular to the flow than in the direction parallel 
to it. This assumption is valid only for unsvept plan forms, which is the 
only case to be considered here. The basic potential equation may there- 
fore be written as 


acL 

AD 

L 

M 

p 
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(l-M^XX + <Pyy + CP ZZ = (r+DM^xCpXX (l) 

where M is the tree-stream Mach number and <P is the perturbation- 
velocity potential normalized by division by the free-stream velocity. 
This equation retains the essential features of both linearized and 
transonic theory and the similarity rules may be obtained from this 
equation and applied at subsonic, transonic, and moderately supersonic 
Mach numbers. 


Similarity of flov about wings of finite span requires the con- 
stancy of two parameters (M 2 -l)A 2 and 


M 2 -! 


(5 


•V 


2/3 


where 


3(x/e)' 

t is an ordinate -amplitude parameter and F 


m 


is the ordi- 


nate distribution function. In general, t may be used to denote 
changes in profile thickness, angle of attack, and camber. The constancy 
of these two parameters implies the existence of an alternative similarity 
parameter 

3 


A 3 t 


3(x/c) 


•F 


(M) 


which is of fundamental importance. Since the ratio of specific heats . 
7 is constant for a fixed medium the expression ( 7 + 1 ) which appears 
in equation (l) has not been retained in the similarity parameters nor 
in the following similarity rules. 

The similarity rule for the pressure coefficient on the surface of 
the wing may be written in the form 




* L T 3 (x/c) r U’bjj p \ r 5 


M 2 -! 


L T 3(x/ c r VcvJ 


r 3 A y> 

1 1/3 X y 1 

_ T 3(x/c) r \c'b/ 



( 2 ) 


3 /x y\ 

The expression t ^ ( x/c ~) ^ '^c , b ' J represents the slope of the airfoil 


surface and can be expressed as 


3(x/c) 


\t) . ( t/c)-_L- \?(%A + igfi.z') - -s.il 
c b/ 3(x/c) |_ \c b/ t/c \c b/ t/c cj 
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where f 


/x y\ , /x y\ 

Ve^by 


are the airfoil thickneas and camber 


distribution functions with the origin at the leading edge of the root 
chord and the parameter h is proportional to the amount of camber. 

For airfoils having the same thickness and camber distribution functions 
the similarity rule for the pressure coefficient may be expressed in the 
more useful form 


C p = (t/c) 5/3 P 


M g -1 ... . 1/3 a _h_ x y 
_(t/c) 2 ^' A( /C) ' t /c' t/c ; c' b 


(b) 


In this equation a and t/c can be simultaneously interchanged, that 
is, t/c may be replaced by a, and a by t/c. This is possible by 
virtue of the presence of the parameter The form of the function P, 

however, may be altered in the process. 


Generalized Force and Moment Coefficients 


From the similarity rule for the pressure coefficient the general- 
ized egressions for lift coefficient and lift-curve slope are obtained 


Cl T I M -1 , , .1/3 a h 

^T 73 ~ L [(t/c) 2 / 3 ' ( /C) 7 t/c* t/cj 


and 


(V=) 1/3 ^ - Lx 


M 2 -l 


M -1 , , . 1/3 a h 

(t/cW*’ A(t/c) ._ / t7c' t7S 


(5) 


( 6 ) 


When the lift coefficient varies linearly with angle of attack the param- 
eter may be omitted from equation (6). 


The generalized expressions for moment coefficient, pitching- 
moment slope, and center of pressure are 


a 


Cm 

2/3 


M 


M 2 -l 


dCm 

dCL 

C.P. 


L<t/e)**’ A(t/c)1/3 ’ V’ 
^(t/c) 2 / 3 ' Alt ^ ' • TP t To] 

pry. 


(7) 

( 8 ) 

(9) 
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The total drag coefficient may be expressed as 


C D = c Df + ^ C Dp^ min + AC D 


( 10 ) 


The friction-drag coefficient CDf cannot be expressed in terms of the 
similarity parameters since the theory neglects all viscous effects. 

The generalized expression for minimum pressure drag is 


i-BsLala =n 

(t/c ) 5 / 3 ^ 


m 2 -i 


L(t/c ) 2 / 3 


■> A(t/c) 1/3 , 

t/c_ 


(H) 


Although the similarity rule implies a generalized drag-due -to-lift 
parameter —^ 73 * a more useful generalized expression is 


(t/e) 1/3 ^ - A Di[(T7cyW A(t/c) 1/3 , ^ (12) 


and a camparison of ACp with Cl 2 may be obtained from 
(t/c )" 1 / 3 * AQj [ (t y c ya73* A(t/c) 1/3 , 


(13) 


In equations (12) and (13) the ratio 


a 

t/c 


t/c_ 

may be neglected if the drag 


due to lift shows a parabolic variation with angle of attack and with 
lift coefficient. 


A comparison of drag due to lift with the limits for drag with full 
leading-edge suction (elliptic spanwise loading) and drag with no leading - 
edge suction may be obtained by use of the expression 


irA 


ACp - Cpct 


or, when written in terms of the generalized lift-curve slope and assum- 
ing a linear variation of Cl with a. 


t(t/ f/u% ]S - (t/=) V3 ^ < = (t/=) 1/3 < 

jtA(t/c) 1 /3 a 2 


'La 


In the preceding generalized expressions the argument 


M a -1 


t/c 


(14) 

vanishes 


for symmetrical profiles, — . . g / 3 vanishes at M = 1, and in the case 

( t /c ) 

of two-dimensional flow A( t/c ) 1 / 3 does not enter into the similarity 
rules. 
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Use of the Similarity Parameters 


the general- 


The generalized coefficients and derivatives can be applied only 
to experimental data for wings having the same thickness and camber dis- 
tribution functions. For constant values of - 5 - and -S— 

t/c t/c 

ized expressions will remain oonstant and the flow patterns will be simi- 
lar if the Mach numbers are chosen so that J [m 2 -1 [a and 

remain constant. This will be possible only if the wings have equal- 
values A(t/c) 1 / 3 , and thus similarity can exist only if the wing geome- 
tries are related in a special manner. The variations of Mach number 
with thickness ratio and with aspect ratio for several constant values of 

J |m 2 -i|a and are shown in figure l(a). The variations 

between aspect ratio and thickness ratio for several constant values of 
A(t/c ) i ^ 3 are shown in figure l(b). 


According to small perturbation theory the thickness ratio must be 
small but no restriction is necessary concerning the magnitude of the 
aspect ratio. A given value of A(t/c ) 1 / 3 must always correspond to a 
small thickness ratio and therefore a systematic survey of airfoil data 
should cover a wide range of A(t/c ) 1 / 3 values as was done for the 
analysis of this paper. (The values of A(t/c ) x / 3 for these airfoils 
are tabulated in table I . ) 


For wings^having unequal values of A(t/c ) 1//3 the use of the two 
parameters -^(t/c ) 1 ^ 3 suggests two different but essen- 

tially equivalent forms of data correlation, one showing the variation 
of the generalized coefficients with A( t/c ) 1/1,3 for constant values of 

M 2 -l 


the speed parameter 


the other showing the variation with 


M 2_i ft/c) 2 ^ 3 ' 

7 — , ~ . for constant values of A(t/c) 1 ' 3 . For convenience, the corre- 
(t/c ) 2 / 3 m2.! 


lation in this paper will be made for constant values of 


and 


(t/c ) 2/3 

the results summarized for constant values of A(t/c) 1 ' 3 . Thus, the use 
of the similarity parameters permits the multiple families of basic data 
curves for various aspect ratios and thickness ratios to.be summarized 
in a presentation involving only one geometric variable, the parameter 
ACt/c) 1 / 3 . 


The actual flow about wings will always violate to some extent the 
assumptions necessary in the derivation of small disturbance theory, 
which assumes that the flow deflection is everywhere small. The theory 
also ignores the existence of a stagnation region at the leading edge 
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and of a boundary layer. A theory for transonic flow must also consider 
the extent and position of the mixed subsonic and supersonic flow regions. 
Although the essential properties of shock waves and of PrandtlrMeyer 
flow can be approximated in terms of the speed parameter, provided the 
flow deflection angle is small (references 1 and 5) , a successful corre- 
lation of experimental data will be possible only if the boundaries of 
the mixed, flow regions for related airfoils possess similitude of size 
and position. 

The usefulness of the similarity rules can be determined only from 
a systematic survey of experimental data. 


Slender -Wing Theory 


A solution in general is not known for the nonlinear transonic 
potential equation 

(l-M 2 )^ + <Pyy + <P ZZ = (r+l)M 2 cp x cp xx (l) 

For two-dimensional flow the partial derivatives ‘Tlx and 9xx are 
known to become very large as the Mach number approaches 1. However, 
for wings of slender plan form, that is, wings of sufficiently low aspect 
ratio, these derivatives may be considered to be of the same relative 
magnitude as the partial derivatives 9yy and ty Z z and the nonlinear 
term (7+l)M S< Px <: Pxx> which is of the second order, can therefore be 
neglected. Thus, for slender wings the linearized equation 

(l-lf^cPxx + qPyy + tPzz = 0 (^5) 

can be used throughout the transonic speed range. 

For vanishingly small aspect ratios, or for moderate aspect ratios 
at near sonic speeds where the coefficient 1-M 2 becomes small independ- 
ently of 'Pxx* the linearized equation reduces to Laplace's equation in 
two dimensions 


y-yj + 'Pzz 0 

for which solutions are well known. Such solutions have been given by 
R. T. Jones (reference 6) for the case of slender, pointed wings and by 
various authors (see, e.g., references 7 and 8) for slender wing-body 
combinations . The slender -wing theory furnishes the interesting result 
that all the lift is carried upstream of the point of maximum span and 
the center of pressure of a rectangular wing of vanishing aspect ratio 
will, at sonic speed, be at the leading edge. The theory, however, is 
unable to predict the effects of profile thickness and is restricted to 
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the lifting case. A solution cannot be obtained for the zero -lift pres- 
sure drag. 

The theory of reference 6 provides the following expressions for 
the lift and drag coefficients: 


Cl = | Aa 

ACp = ^ Aa 2 (full leading-edge suction) 

(16) 
(IT) 

will be used in the data correlation to obtain a comparison between 
theory and experiment for the wings, of low aspect ratio at the sonic 
Mach number. 


The alternative expressions 

( t / = ) vs 3a, „ | A(t/c ) 1/3 

da 2 

(t/c) 1 ^^^ = ^ A(t/c) ly/3 (full leading-edge suction) 


DATA CORRELATION 


Before beginning the data correlation it might be well to point 
out that the experimental data were obtained by mounting semispan wing 
models in the high-velocity flow field of the Ames l6-foot high-speed 
wind-tunnel tr ans onic bump. .The streamlines of the bump flow field are 
slightly curved with Mach number gradients in the plane of the wing 
model. Typical Mach number contours are presented in reference 4. 

These gradients are most pronounced at the higher Mach numbers and for 
the larger aspect ratios . The effects of the nonuniformity of the flow 
field are unknown but a certain rounding off of any sharp breaks in force 
and moment coefficient variation with Mach number can be expected. 

Minimum Pressure Drag 


The drag-similarity rule cannot be applied directly to a correla- 
tion of minimum drag data since small perturbation theory ignores the 
existence of friction drag. The friction drag coefficient is believed 
to change little with Maph number in the transonic range and the corre- 
lation can be applied to the minimum pressure drag coefficient. The 
minimum pressure drag has been calculated by subtracting a constant, 
friction drag, assumed equal to the minimum drag at 0.7 Mach number, 
from the minimum drag, that is, " 


(^Dp)min ^Qmin “ ^Dmin^M=o.7 


(18) 
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The basic data curves for the variation of minimum drag coefficient 
with Mach number are presented in figure 2. The variation of the general- 


ized coefficient 


(t/c) 5 / 3 


with the speed parameter 


M 2 -! 

(t/c)2/3 


is shown 


in figure 3 as an introductive step to the completely generalized data 
correlation of figures 1(- and 5- (The symbols used in fig. 3 represent 
the actual test data points and were included to show the manner in 
which the data have been faired. ) 


The data curves for the wings of aspect ratio 6 in figure 3 niay tie 
closely represented by a common curve, indicating that at this aspect 

(Cp-) 

m n may be considered to depend only on the 


ratio the values of 
M 2 -l 


value of 


(t/c) 2 / 3 


(t/c)S/ 3 

and that the flow is essentially two-dimensional. 


The aspect -ratio -6 wings also exhibit the negative variation of the 
force coefficient with Mach number which is characteristic of two- 
dimensional flows at sonic speed. This variation is a consequence of 
the relative variations of local and free-stream dynamic pressures. 

The local Mach numbers are effectively frozen at near -sonic values as 
the free-stream Mach number increases through the transonic speed range. 
(See references 9 and 10.) It is interesting to note that the experi- 
mental values of ( ■ dCl) ) for the wings of aspect ratio 6, figure 2, 

V dMyta =1 

agree qualitatively with the values implied by the following relation- 
ship (an exact theoretical result applicable to symmetrical profiles 
of any shape, reference 10) : 



2 

7+1 


(°»W 


(19) 


This agreement occurred in spite of the fact that the test conditions 
are not ideal and do not agree with the concept of an infinite and 
uniform flow field assumed in the theoretical reasoning. 


The lower aspect ratios show increasing effects of three-dimensional 
flow and the smaller thickness ratios have progressively lower general- 
ized drag coefficients with the exception of the thinnest (2-percent- 
thick) wing models. These thinnest airfoils have unusually large values 


of 


( ^Dp ^min 
(t/c) 5/3* 


which are believed to be largely the result of the boundary 


layer creating an effective thickness considerably larger than the 
actual profile thickness. The generalized drag coefficients are then 
magnified for the thin airfoils by the 5/3 powers of the ratios of 
effective thickness to profile thickness. 
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The completely generalized correlation of minimum pressure drag is 
presented in figures 4 and 5. (The symbols used in these figures and in 
the remaining figures of the report represent values taken from faired 
curves of preliminary cross plots. A different symbol has been used to 
represent the data for each particula r wing t ested.) The transonic 

correlation for constant values of J | 1-M a j A is presented in figure If 
merely to illustrate the inability of presenting data at the sonic 

velocity using simultaneously the parameters *J |m 2 -1|a and - — — 2 / 3 • 

M s -1 ' ' c 

The correlation for constant values of Presented in 


figure 5 . The correlation is best at the sonic velocity and the curve 
for this particularly interesting Mach number is repeated in figure 6. 
Although the data for 2-percent thickness have been omitted from the 
preceding correlation because the extraordinarily large values interfere 
with the adjacent correlation curves, these data have been included in 
figure 6 to show the pronounced effect of the boundary layer for these 
thinnest profiles. 


At the sonic speed the minimum pressure drag is seen to vary lin- 
early with aspect ratio and with the second power of the thickness ratio 
for values of A(t/c) 1//3 less than about 1 as can be expressed by 




M 


A(t/c) 1/3 C 1 


( 20 ) 


For values of A(t/c) 1 / 3 greater than about 1 the generalized coeffi- 


cient 


(Cd pU 


— — y s yS approaches rapidly and asymptotically toward a constant 
(t/c ) 5 / 3 

value for which the minimum pressure drag varies with the 5/3 power of 
the thickness ratio in accordance with the drag similarity rule for 
sonic, two-dimensional flow. The extrapolated two-dimensional-flow 
value for M=l, (CDp) min = 3*55 (t/c) 5 / 3 , was obtained by plotting 

against the inverse parameter l/A(t/e) 1//3 . The calculated theoretical 
pressure drag for a double-wedge profile (reference ll) is somewhat 
higher . 


The correlation. of minimum pressure drag is summarized in figure 7 
by cross plotting from the faired curves of figure 5- The predicted 
critical Mach numbers 1 for MCA ^30XX profiles are given in reference 13 

and may be approximated by 0 /o = -1.95* The compressibility drag 


(t/c) S/3 


T In reference 12 Kaplan has shown that the section critical Mach number, 
according to first-order linearized theory, is related to the thick- 
ness ratio by a constant value of the speed parameter, the constant 
depending on the particular profile shape. 
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rise becomes noticeable at slightly higher Mach numbers, that is, for 
M g -1 _ , o n 
(t/c ) e / 3 8 


The fundamental importance of the similarity parameter A(t/c ) 1 / 3 
is clearly evident from the summarization of minimum pressure drag as 
presented in figure 7 . The curves of figure 1 and the values of table II 
may be used- with figure 7 for convenience in calculating the minimum 
pressure drag for given values of aspect ratio, thickness ratio, and 
Mach number . 


Lift 


The variation of the lift coefficient with angle of attaek is 
essentially linear at moderate angles of attack throughout . the Mach 
number range for aspect ratios greater than 1.5. The lower aspect 
ratios, however, show an increasingly nonlinear variation of lift with 
angle of attack and approach the theoretical sin? a variation for 
vanishing aspect ratio. For convenience the lift analysis will be 
restricted to a consideration of the lift-curve slope evaluated at zero- 
lift coefficient, which provides a close approximation for lift charac- 
teristics at the moderate angles of attack for which the s imi larity 
rules can be expected to hold. 

The variation of lift-curve slope with Mach number is shown in 
figure 8 and Is compared with the theoretical lift-curve slopes calcu- 
lated by applying the three-dimensional Prandtl-Glauert transformation 
to the Weissinger lifting-line theory of reference H. The agreement 
between theory and experiment for aspect ratios greater than about 3 is 
satisfactory only in the subcritical Mach number range. 

Above the critical Mach number an abrupt decrease in lift- curve 
slope occurs for some of the wings of larger aspect ratio and thickness 
ratio and is believed to be the result of the formation of strong 
velocity discontinuities and flow separation at the airfoil surfaces . 

This "bucket" type variation in the lift-curve slope is a phenomenon 
apparently dependent on a combination of thickness and aspect-ratio 
effect since the smaller thickness ratios and smaller aspect ratios 
show no such irregularities in the lift-curve variation with Mach number. 
Indeed, It will be shown in the following data correlation that this 
erratic variation occurs only for the wings having values of A( t/c ) 1 / 3 
greater than about 1 . 6 . 

The variation of the generalized lift-curve-Blope parameter 
(t/c ) 1/3 with the speed parameter Is shown In 


djOT KIJlHJEEIA L 
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figure 9 as an introductive step to the completely generalized data 
correlation of figure 10. A comparison of figures 8 and 9 shows that 
the curves for aspect ratio 6 have converged toward a single variation 
at the transonic speeds, indicating that the flow is essentially two- 
dimensional for the higher aspect ratios. „For these wings the lift- 
curve slope shows the characteristic negative variation with Mach number 
at near-sonic speeds. 


of 


The correlation of (t/c) 
M 2 -l 


1/3 


/dCL\ 

\ da / a-o 


for several constant values 


is presented in figure 10. The correlation and comparison 


(t/c) 2/3 

with theory are good for the large negative values of the speed param- 
eter (i.e-. , low Mach numbers) and a reasonably good correlation is 
indicated at near-sonic speeds. At supercritical Mach numbers and 
values of A( t/c ) 1 / 3 greater than 1.6 the poor correlation suggests a 
separated flow for which the concepts of small perturbation theory 
would be violated. The value A(t/c) 1 / 3 = 1.6 appears to be the limit- 
ing value for airfoils which do not exhibit noticeable irregularities in 
lift -curve -slope variation with Mach number. 


At M=l, and for small values of A(t/c) x / 3 , the experimental lift 
agrees well with the lifting-line and slender -wing theories. Hence, 
within the indicated limits, the lift may be approximated by' the follow- 
ing equation: 



M = 1 
A( t/c) 1/3 < 1 


( 21 ) 


For increasingly greater values of A( t/c) 1 ^ 3 a rapid and apparently 

asymptotic approach of (t/o)Vs(^M a _^ to a constant value is 

indicated. The theoretical result of Guderley and Yoshihara 2 for a 
double-wedge profile in sonic, two-dimensional flow is included in 
figure 10. 


The results of the data correlation for lift-curve slope are summa- 
rized in figure 11 for those values of A(t/c) 1 / 3 and ^/ e ' ^g/ a "where 
a good correlation was indicated. 


2 This theoretical value of the lift -curve slope was calculated in 

Air -Force Technical Report 6683, Unsymmetric Flow Patterns at Mach 
Number One, by Gottfried Guderley and Hideo Yoshihara. 
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Moment 


The variation of the degree of static longitudinal stability with 
Mach number is presented in figure 12. The pitching -moment-curve Blope 

dC m /dCL was measured for ratios of 0, 1, and 2 radians. (The 

t/c 

values of a in degrees are tabulated in table III.) 

For the larger aspect ratios varies smoothly with Mach number 

dCL 

in the low subsonic and in the supersonic regimes, although at different 
levels, but behaves erratically at high subsonic speeds. This erratic 
behavior is probably due to the same causes as the irregularities of the 
lift-curve slope at the corresponding Mach numbers. At Mach numbers up 
to the critical the effects of compressibility are relatively small as 
is predicted by linearized theory. At transonic speeds a large change 
in the center of pressure occurs for the wings of large aspect ratio as 
the aerodynamic center moves from the vicinity of the 25 -percent-chord 
point of subsonic speeds, towards the 40-percent -chord point. Only for 
very low values of the aspect ratio is this undesirable change in center 
of pressure substantially decreased. 


The correlation of pitching -moment-curve slope 

M 2 -! 


dCni 

dC L 


for various 


values of the speed parameter is presented in figure 13 . 


Although considerable scatter of data is evident the curves have been 
faired favoring the data for the thickness ratio of 4 percent which do 

show a good correlation. The values of — ^ for = 0 represent 

dCL t/c 

the position of the center of pressure and for all values of the speed 
parameter the center of pressure is shown to move progressively toward 
the leading edge as the aspect ratio becomes small. At transonic speeds 
the stability derivative varies almost linearly with A(tZc) 1 / 3 for 
values of A(t/c) 1 / 3 less than about 1 with the center of pressure 
located at the leading edge for vanishing aspect ratio at zero angle of 
attack. For values of ACt/c) 1 / 3 greater than about 1 the stability 
derivatives may be considered constant and Independent of both aspect 
ratio and thickness ratio. 


The pitching-moment-curve slope correlation is summarized in fig- 
ure 14. When using these curves to estimate — - for particular values ' 

CU,L 

of aspect ratio and thickness ratio the lift coefficient corresponding 

a, 

to a given value of can be approximated by use of the identity 


Cl = (t/c) 2 / 3 (t/c) 1 / 3 ^ 


t/c 
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In this equation the lift variation with angle of attack has been 
assumed to be linear. Numerical values for ( t/c) 1 / 3 Cy i can be 

obtained from the s ummar y curves of figure 11. . 

Drag Due to Lift 

Equation (12) will be used for the correlation of the drag-due-to- 


lift coefficient 
stant values of 


ACj-j. The correlation will be made for several con- 
in radians), since it is not always known 


beforehand if it is possible to represent closely the drag due to lift 
by a parabolic variation with the angle of .attack. 


The variation of ( t/c vith the speed parameter 


M a -1 


a 


(t/=) 2/3 

for the various thickness ratios and aspect ratios is shown in figure 15 . 
The approximate limits for drag with full leading-edge suction (elliptic 
spanwise loading) and drag with no leading-edge suction as given by 
equation ( Ilf) have been evaluated using the summary curves of figure 11 
for the lift-curve slope and are presented in figure 15 to show the 
degree of leading-edge suction. At transonic speeds the drag force for 
the larger aspect ratios and thickness ratios is actually somewhat 
higher than the value corresponding to a resultant force perpendicular 
to the plane of the wing, suggesting that some increase lp separation 
and viscous effects occurs with increasing angle of attack. 

The correlation of (t/c ^3 for several constant values of 

m 2_i _ a a 

(t/c) g/s t/c P resen ^ e< ^ * n figure. 16. The correlation curves 

for constant — Si- are presented in the left-hand side and summarized 

Vc 

in the right-hand side of figure 16. A poor correlation is Indicated 
for large Values of A(t/c) 1 ' 3 at the large negative values of the speed 
parameter where the degree of leading-edge suction .apparently is chang- 
ing rapidly with Mach number. When the resultant force becomes normal 
to the chord line the drag -due -to -lift correlation is connected inti- 
mately with the lift correlation. For transonic values of the speed 

parameter a reasonably good correlation of ( t/c ) 1 ^ 3 is indicated 

a 

and the various curves may be closely represented by a common 

curve as they should for a parabolic variation of ACj) with a. 


The sonic data are presented in figure 16(c) and are compared with 

Ct® 

the low-aspect-ratio theory for both ACp = -**- and ACj) = Cm (where 

jtA 

the lift coefficient is given by Cl = 75 Aa) . These s ummar y curves show 
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that the drag due to lift, for low-aspect ratios, is given approximately 
by 


ACp = | Aa z ; 


M = 1 

0.5<A(t/c) 1 / 3 <1.3 


( 22 ) 


Only for very low values of A(t/c) 1 ^ 3 does the drag due to lift tend 
toward the formal result of low-aspect-ratio theory, ACp = - Aa 2 . 


The drag-due- to -lift correlation is summarized in figure 17- Since 
the correlation curves were found to be nearly independent of 

especially in the transonic speed range and for the moderate angles of 

cc. 

attack used, the presentation is made only for the ratio ^ - = 1. 

A comparison of the induced drag for constant lift coefficient may 
be obtained from a consideration of the generalized parameter 

(t/c) 1/ ' S ^2. The variation of this parameter with — M 2 -!- — f or con- 

CL 2 , 1/a (t/c) 2 /a 

stant values of A(t/c) /s is shown in figure 18 as derived from the 
lift and drag summary curves of figures 11 and 17* A consistently 
decreasing drag due to lift, for constant lift coefficient, is indicated 
for increasing aspect ratio. It should be remembered, however, that the 
correlation of data does not apply at the high subsonic Mach numbers for 
values of ACt/c) 1 / 3 greater than 1.6. For these wings, of large 
aspect ratio and large thickness ratio, a reversed and erratic trend of 
large induced drag results, probably from the large effects of shock- 
induced flow separation. 


ADDITIONAL CONSIDERATIONS 

Maximum Lift-Drag Ratio and Optimum Lift Coefficient 


The variation of the maximum lift-drag ratio and the corresponding 
opt imum lift coefficient with Mach number is shown in figures 19 and 20 . 


The maximum lift-drag ratio is expressed by the familiar formula 

(0 


max 



[0D f + (CBpW 


(23) 


provided the variation of ACp with Cp is parabolic. The value of the 
the friction coefficient CPf for the conditions under consideration 
may be taken to be 0.006. From the equ ivalen t expression 


Mb’ mkltfl'TM . 
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L(V^ F7cf7 5 'J (t/c) 


ns 


v-l/3 


AC n 

cl 2 


(24) 


it is apparent that the inclusion of the friction drag does not permit 

a generalized expression for f —) in terms of the similarity param- 

\D /max 

eters. A correlation of data must cons ider # individually, each separate 
value of the thickness ratio. Although no generalized correlation 
appears possible a few remarks will be made concerning the effects of 
aspect ratio and thickness .ratio at transonic speeds. 


At low speeds and large Reynolds numbers, the maximum lift-drag 
ratio shows a dependence principally on the value of the aspect ratio. 
Above the critical Mach number, where the minimum pressure drag becomes 
large, the lift-drag ratios are found to be essentially independent of 
the aspect ratio but indicate a pronounced dependence on the. thickness 
ratio. This independence of aspect ratio at transonic speeds is illus- 
trated by the sonic data in figure 21 where the variation with A(tfc) 1 / 3 
of the two expressions 


' C Df - 


L (t/c) 5/3 j 



ACp 


(t/c)' 


■1/3 

and 

f ( CD p)min 


- 


L (t/c)s/aj 



■^£ (t/c)’ 1 ' 3 ] 


used in equation ( 23 ) for f g ie shown for several thickness 

ratios. The two variations have a compensating effect that leads to a 
maximum lift-drag ratio essentially independent of aspect ratio. 


The optimum lift coefficient is represented by the following for- 
mula when the drag polar s are parabolic: 

( c Dp)inin 

ACp V ^ 

Cl 2 

The basic data curves of figure 20 show that the value of CLop^ at 
near-sonic speeds decreases with both decreasing aspect ratio and 
decreasing thickness ratio. 



The linear variation of the force coefficients with aspect ratio 
at the sonic speed is characteristic of the. wings having values of 
A(t/c) x / 3 less than about 1. The following empirical formulas were 
obtained by substituting equations (20), (21), and (22) into equa- 
tions ( 23 ) and ( 25 ) 
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Hie friction drag coefficient may be closely approximated by 0.006 and 
the calculated values from, these formulas agree closely -with the sonic 
experimental data for the aspect-ratio range of 3 to 1. 


Drag-Divergence Mach Humber 


The drag-divergence Mach number Mpp may be defined as the Mach 
number at which a definite and abrupt increase in drag coefficient 
occurs and is usually chosen as the Mach number for which the rate of 
change of drag coefficient with Mach number reaches some arbitrary 
value, say 0.1. A new definition for drag-divergence Mach n umb er follows 
from a consideration of the similarity rule for the drag coefficient 


CD = (t/c ) 5/3 D 


r m 2 -i 

_(t/c)2/3’ 


A(t/c) 1/S , 


q h, ~| 
t/c* t/c J 


Differentiating with respect to the Mach number gives 


- (t/c)M Di 


(t/c) 2 / 3 ’ A ( t / c ) 1/3 j 


a h 
t/c’ t/c 


The drag-divergence Mach number may now be defined as the Mach n umb er 
for which 


dCp 

dM 

(t/c) M 


= constant 


and this implies the functional relationship 


l-Mpp _ 

(t/0)*' 8 ' 


A(t/o)1/3 > t7?’ t7S 


( 27 ) 


( 28 ) 


For two-dimensional flows it follows immediately that the relationship 
between drag-divergence Mach number and thickness ratio for symmetrical 

l-Mp© 2 


wings at ?ero angle of attack is given by a constant value of 


(t/c) 2 / 3 
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The_ experimental drag-divergence Mach numbers for zero lift coeffi- 
<x \ 

r~r = 0 are correlated in figure 22. For convenience the 

b / G / acD/aM 

evaluation has bben made for = -*-■ According to figure 22 the 

two -dimensional -flow value for Mdd is given approximately by 

= 1*75 as is indicated by the asymptotic nature of the curve. 


CONCLUDING REMARKS 


The similarity rules have been used to correlate the experimental 
data for a series of 22 rectangular, symmetrical wings having 
NACA 63AOXX sections, aspect ratios from l/2 to 6, and thicknesses frcan 

2 to 10 percent. The data were obtained by use of the transonic bump 
technique for a Mach number range of 0 .40 to 1.10 and Reynolds number 
range of 1.25 to 2.05 million. 

The results of the correlation have Bhown that, with the exception 
of wings having large values of A(t/c) 1//3 at high subsonic Mach numbers 
where an erratic variation of the force and moment coefficients with 
Mach number was indicated, it is possible to correlate experimental data 
throughout the Mach number range using the transonic similarity param- 
eters. The use of the generalized coefficients has permitted the pres- 
entation of experimental data for a wide range of aspect ratios and 
thickness ratios by -a unified method throughout the Mach number range, 
and the form of presentation used has permitted a direct comparison of 
the data with the known results of theory. 

At the sonic Mach number a linear variation of the force and moment 
coefficients with aspect ratio was found to be a universal property for 
wings having values of A( t/c) 1 / 3 less than about 1. For increasing 
values of A(t/c) a / 3 greater than 1 the generalized coefficients at 
the sonic speed show a rapid "and asymptotic approach to constant values. 
Indicating that a transition from three-dimensional-flow characteristics 
(where the force and moment coefficients vary linearily with aspect 
ratio) to two -dimensional -flow characteristics (where the force and 
moment coefficients are essentially independent of the aspect ratio) 
occurs near the particular value of A(t/c ) 2 / 3 equal to 1, 

The data correlation was summarized in presentations involving only 
one geometric variable A(t/c ) 1 / 3 . The summary curves may be used as 
design charts for estimating the transonic characteristics of rectangular 
wings provided the airfoil profile does not differ greatly from the 
NACA 63 A series. Although a correlation of experimental data was not 
possible for the maximum lift-drag ratio and the optimum lift coefficient, 


mm 



NACA RM A51L17b 


the summary curves for minimum pressure drag and drag due to lift may 
be used to estimate these values. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics , 
Moffett Field, Calif. 
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TABLE I.- VALUES OF THE GEOMETRIC PARAMETER 
ACt/c) 1 / 3 FOR THE WING MODELS TESTED 


A 

t/c 

A(t/0)V3 


A 

t/c 


6 

0.10 

2.784 


2 

o.o4 

0.684 

6 

.08 

2.586. 


2 

.02 

.542 

6 

.06 

2.352 


1.5 

.04 

.513 

4 

.10 

1.856 


1-5 

.02 

.407 

4 

.08 

1.724 


1 

.10 

.464 

4 

.06 

1.568 


1 

.08 

.431 

4 

.o4 

1.368 


1 

.06 

.392 

3 

.04 

1.026 


1 

.04 

, 3^2 

2 

.10 

.928 


1 

.02 

.271 

2 

.08 

.862 


• 5 

.04 

.171 

2 

.06 

.784 


• 5 

.02 

.136 


<*sote§^£l^ 


4R 
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'CABLE II.- NUMERICAL VALUES OF (t/c) 1/3 , 
(t/c) 2/3 , AND (t/c) 5/s 


t/c 

(t/c) 1/a 

(t/c) 2/3 

(t/c) 5 / 2 

0.12 

0.493 

0.243 

0.0293 

.11 

.479 . 

.230 

.0254 

.10 

.464 

.216 

.0217 

.09 

.448 

.201 

.0181 

.08 

.431 

.186 

.0148 

. .07 

.412 

.170 

.0118 

.06 . 

.392 

.154 

.0092 

.05 

.368 

.136 

.0068 

.04 

.342 

.117 

.00468 

.03 

• 311 

.097 

.00292 

.02 

.271 

.073 

.00147 
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TABLE III.- VALUES OF THE ANGLE OF ATTACK (IN DEGREES) 

FOR VARIOUS VALUES OF THE RATIO ~~ (IN RADIANS) 

t/c 


t/c 

■Angle of attack a, degrees 

a , 0.5 radian 
t/c ' 

1.0 radian 

1.5 radians 

2.0 radians 

0.10 

2.87 

5-73 

8.60 

- - 

.08 

2.29 

4.58 

6.87 

- - 

.06 

1.72 

3-44 

5.16 

6.87 

.04 

- - 

2.29 

3-44 

4.58 

.02 

- - 

1 

1 

1.72 

2.29 
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Aspect ratio, A Thickness ratio, t/c 

(a) The functional variation of aspect ratio and thickness ratio 

with Mach number. 

Figure i.~ Requirements for similitude of transonic flow. 
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Figure 3.- The variation of 


(Qop)mln 
t t/c) V* 


with 


M*-i 

(t/c)*' 3 


for the various aspect ratios and 


thickness ratios. 
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Figure 5.- The correlation of the generalized minimum pressure-drag coefficient 
for constant values of the speed parameter. 
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Figure /5.- The variation of the generalized drag-due-to-iift para- 
meter with the speed parameter for several 
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ratios of a/(t/c). a in radians. 
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Figure 20.- 


The variation 


of lift coefficient for maximum tift-drag rotio with 


Mach number . 
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